Phase change in an opinion-dynamics model with separation of time scales 
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We define an opinion formation model of agents in a Id ring, where the opinion of an agent evolves 
due to its interactions with close neighbors and due to its either positive or negative attitude toward 
the overall mood of all the other agents. While the dynamics of the agent's opinion is described 
with an appropriate differential equation, from time to time pairs of agents are allowed to change 
their locations to improve the homogeneity of opinion (or comfort feeling) with respect to their short 
range environment. In this way the time scale of transaction dynamics and that of environment 
update are well separated and controlled by a single parameter. By varying this parameter we 
discovered a phase change in the number of undecided individuals. This phenomenon arises from 
the fact that too frequent location exchanges among agents result in frustration in their opinion 
formation. Our mean field analysis supports this picture. 

PACS numbers: 89.75.Fb, 87.23.Ge, 64.60.aq 



I. INTRODUCTION 

How are opinions formed? In sociology this is one of 
the basic questions, but it is also highly relevant for poli- 
tics, innovation spreading, decision making, and the gen- 
eral well feeling of people [TH3] • This complex process de- 
pends on various factors or components like confidence, 
attitudes, communities or media effects [4j. Recently, 
much effort has been invested in modeling different as- 
pects of opinion dynamics and these models are in many 
ways related to those of physics [6]. Unfortunately, 
the empirical observations are rather sparse. Therefore, 
the usual strategy is to concentrate on some particular 
features by making plausible assumptions for a model, 
and comparing its results with expectations. Here we 
will follow this line of study. 

Our starting point is that the comfort feeling of an 
individual depends on his/her embedding in the soci- 
ety. We get friends mostly with people who are simi- 
lar to us, share our opinions, tastes etc. In sociology 
this is called homophily and is known to be the major 
governing principle in friendship formation pQ. In terms 
of physics, this corresponds to ferromagnetic interactions 
[7j. In the language of opinion dynamics this means that: 
a) The opinion of an individual gets adjusted to that of 
his/her friendship neighborhood; b) An individual seeks 
the neighborhood of alike others. Here a) has been the 
basis of many opinion-dynamics models, both discrete 
and continuous [5H14). On the other hand b) has been 
investigated in the framework of coevolving networks |15l - 
I21j . where the connections between individuals are not 
there forever but can be changed in parallel with the 
evolution of the opinions in order to increase the level of 
satisfaction in the system [52H25]. 

Recently we have introduced a coevolving network 
model [281 129] , where not only short range ferromagnetic 
interactions but also long range interactions were taken 
into account. This corresponds to the fact that, although 



our opinion is strongly influenced by our close friends, we 
are not independent of the general mood of the society. 
However, the impact of the society as a whole does not 
have to be ferromagnetic. As known from sociology again 
PQ , all individuals have two kinds of driving forces with 
respect to the society: We want to be similar to the aver- 
age around us to use the society's protecting power and, 
at the same time, we want to be different to be distin- 
guished as individuals. For every individual these con- 
flicting components are present in different proportions, 
resulting in either net positive or net negative attitude 
with respect to the overall opinion of other individuals. 
This effect was taken into account [28] by an attitude 
parameter a considered fixed or quenched to each indi- 
vidual. Since the attitude parameter can have positive 
or negative sign, it constitutes a source of frustration |30] 
in the system. 

In [2H1 HH] we also introduced a separation of time 
scales for different opinion formation mechanisms. While 
communications go on all the time leading to a quasi- 
continuous adjustment of the individuals' opinions, it 
takes more effort to make new friends than to quit with 
old ones. Therefore, we introduced a measure of time 
separation, which characterizes this difference by allow- 
ing for changes in the network neighborhoods after g time 
steps of the difference equation governing the opinion up- 
date. We found interesting effects as a function of g and 
the attitude parameter a: For small values of g, where 
the rewiring process is very rapid and only two commu- 
nities eventually develop, the attitude parameter plays a 
minor role and the a distributions in them were found 
to be broad and similar. However, for the intermediate 
values of g the smaller communities have a rather narrow 
distribution with mostly negative a values, while the dis- 
tributions for larger communities are broad and shifted 
toward positive a values. Naturally, the agents with neg- 
ative a's do not feel comfortable in a large homogeneous 
community, thus they tend to build smaller ones. 
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The aim of the present paper is to understand bet- 
ter the role of the attitude parameter and the separation 
of time scales in the coevolution of opinion and network 
structure of the underlying system. In order to do so, 
we define a model on a ring, and keep this topology pre- 
served. Therefore, instead of rewiring the network we 
allow for location exchanges between agents by carry- 
ing their individual opinions and attitudes. This corre- 
sponds to a situation where the agent looks for a better 
environment to live in, and is reminiscent of Schelling's 
checkerboard model for residential segregation, where the 
relocation of agents with a mild preference for having a 
few alike neighbors in a static lattice can lead to fully 
segregated outcomes [3TH37] . The decision whether such 
an exchange is made is assumed dependent only on the 
short range interactions. However, in the opinion forma- 
tion the attitude toward the social mood plays also a role. 
As the possibilities for finding new environments are lim- 
ited, an amount of frustration will remain in the system 
for not too large values of the time separation parameter 
g. Interestingly, we see as a function of g a rather sharp, 
phase transition-like change to a state without frustra- 
tion as the individuals get enough time to form a firm 
opinion. 

This paper is organized as follows. In Section [il] we in- 
troduce the model in detail. In Scction [nT| wc present the 
numerical results. In Section HVl a mean field calculation 
is presented, giving account for the variation in the num- 
ber of undecided agents. Finally we draw conclusions. 



II. MODEL 

As in [2"5l |2"5] , we study the dynamics of opinion for- 
mation in a network with a fixed number of individu- 
als or agents (JV) to whom a simple question is posed. 
For the network connectivity between agents, we here 
assume a Id ring topology instead of a more complex 
network topology we studied earlier. A state variable 
Xi € [— Xum, xum] (for fixed xu m > 0) is associated with 
each individual i, which measures the agent's instanta- 
neous inclination concerning the question at hand, while 
the network links represent the presence of discussions 
between agents related to this question. The time scale 
for discussions or exchange of information between indi- 
viduals ("transactions") is dt, while the time scale for a 
generalized change of connections in the network ( "gener- 
ation" ) is T. These two quantities are related by T — gdt, 
where the parameter g defines the number of transactions 
per generation. 

The dynamics of the agent's state variable Xi can be 
written as 



dxj 
dt 



= Is{{x J } s )x i + fi({xj}i)ai 



(1) 



where the random parameter <Xj G [— otu m , &Um] (for fixed 
a Um > 0) accounts for the agent's own attitude towards 
overall or public opinion. The short range interaction 



term f s ({xj} s )xi describes the direct influence over i of 
the subset of 'close' agents {xj} s , while the long range 
interaction term fi({xj}i)aii measures the indirect effect 
of the subset of 'far' agents {xj}i modulated by the at- 
titude of i. The system consist of a ring (a chain with 
periodic boundary conditions) where the short range in- 
teractions take place over the first m neighbors of each 
agent, so the number of short range connections is 2m. 
The long range interaction takes into account the average 
of opinion over the rest of agents in the network, that is, 

i _ j m 

fs{{xj} s )xi = (x)- m) sgn(xi)xi = -^-^x i± t, (2) 



[N/2] 



fl({xj}i)ai = (x)\ 



"/ = ~ : ; Xi±t, 



N - 1 - 2m 



t=m+\ 



(3) 



where sgn(xi) denotes the sign of X{. Observe that m < 
(N — l)/2. Once the opinion of an agent reaches any 
of the limit values ±xu m , it stays fixed for the rest of 
the dynamics and the agent is said to be decided. This is 
because we attempt to describe a state of total conviction 
that is unlikely to change anymore, like in a balloting 
process. 

The dynamical evolution of the systems obeys Eq. [T] 
for g time steps, when the agents are allowed to exchange 
places in the ring in order to help them reaching a defi- 
nite decision (\xi\ = xu m ). This is done according to the 
following rules: One chooses iV 2 pairs of agents at ran- 
dom, and picks up the pairs of agents with both of them 
being not decided (|xj|, \xj\ < xu m ). For these pairs one 
calculates a measure of the distance between the agents' 
opinions 
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(4) 



and compares it with the same quantity if one exchanges 
i and j, namely 



<l 
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If p[" 1 ^ > one exchanges places. In the above 

( - m ' ) means the average over the m left and 



formulas 

right neighbors of site i, as in Eq. [2] This procedure is 
repeated every g time steps, until one is not able to find 
favorable changes or all the agents have reached either 
one of the two limit opinions, since the exchange process 
only deals with pairs of undecided agents. Observe that 
these rules tend to increase opinion homogeneity in the 
system, which is reminiscent of the homophily principle 
mentioned in the previous section. A descriptive diagram 
of the system and the exchange process is shown in Fig. [T] 
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FIG. 1: (Color online) Diagram showing the exchange process 
used in the model (for m = 1). The randomly chosen unde- 
cided nodes i and j have < p^V, thus they are exchanged. 
Observe that the grey scale representing the opinion variable 
is more uniform after exchange. 



III. NUMERICAL RESULTS 

We solve the model system by numerical simulations. 
For that the system is initialized with values of the state 
variable Xi(0) in the interval [—xn m , xu rn ] drawn ran- 
domly from the Gaussian distribution with zero mean 
and unit standard deviation, cut off at ±xu m and Xu m = 
1. Likewise the attitude parameter on for each agent was 
chosen randomly from a uniform distribution in the inter- 
val [— a.umi a Um] with aum = 1 and kept fixed through- 
out the whole simulations. 

The simulations have been carried out using a two- 
step process: first we solve the dynamical equations of 
opinion for all undecided agents by using a simple Eu- 
ler numerical integration with time step dt = ICP 4 , and 
then we perform the exchange process every g time steps 
according to the above described rules. We keep track of 
the progress of the dynamics with two counters, namely 
the number of undecided agents (n un d) and the number 
of pairs that have been exchanged (n exc h) after every 
exchange process. As the agents reach the definite opin- 
ions, the counter n un d will decrease from its initial value 
N to some number close to zero. The asymptotic sta- 
tionary value of n un d is considered as the final number of 
undecided agents. Since the exchanged pairs have to be 
undecided, the counter n exc h usually stays around or be- 
low n^ nd . The exchange process is realized sequentially 
and randomly such that agents can be chosen more than 
once in the same generation. However, the probability of 
such event decays fast with N to be very rear to have an 
effect on the results. 

In the simulations the dynamics is let to run until the 
exchange of agent locations takes place very rarely. The 
relaxation time for this is exceedingly large, and com- 
paring the results of calculations with a large number of 
iterations we found that after 10 7 transactions, the aver- 
aged results over 100 realizations differ by less that 0.1%. 
Therefore, in all the calculations presented here we have 
used these numbers. Moreover, since it turned out that 
some results depend strongly on the size of the ring for 



small values of N, we chose to do most of the calculations 
on a ring of N — 5000 and for the case m = 1, i.e. the 
short range being limited to nearest neighbors. Finally, 
the exchange process rules can be relaxed so that any 
agent (decided or undecided) can be moved, yet since the 
exchange of decided agents is so seldom as to have any 
noticeable effect on the averaged results, we chose not 
to do so in order to increase the speed of the algorithm, 
especially for late stages of the dynamics. 
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FIG. 2: (Color online) (a) Average number of undecided 
agents as a function of the number of transactions per gen- 
eration g. The purple squares are the numerical results after 
10 7 transactions and the red circles are the agents who will 
never get a decision, according to a linear analysis. The corre- 
sponding mean field predictions are shown as continuous lines, 
(b) The average fraction of undecided agents with negative a 
(purple squares) and its corresponding mean field prediction 
(continuous line), (c) Average cluster size as a function of g. 
Observe the plateau around g c « 1.7 X 10 3 . 

In Fig. [2ja) we show the average number of undecided 
agents {{n un d)) in the ring as a function of the parameter 
g. It is clearly seen that there is a quite sharp minimum 
at g c R3 1.7 x 10 3 , a critical value of g that can be pre- 



dicted by using the mean field analysis, see SectionpA 7 ] In 
the figure we also show by red circles the expected value 
of (n un d) f° r t — ¥ oo as obtained from a linear analy- 
sis, and by continuous lines the corresponding mean field 
predictions, explained below in more detail. The criti- 
cal value g c signals a change of phase in the system: for 
g > g c all the agents get decided in the limit of infinite 
time, while for g < g c a finite fraction of the network re- 
mains undecided for arbitrarily long times. We identify 
the former phase as a state of maximum relaxation, and 
the latter as a frustrated state where many agents cannot 
reach the limit opinions. 
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In Fig. [2jb) it is seen that the average value of the 
fraction of undecided agents with negative a ((fa)) a ls° 
shows a sharp change of behavior, predicted by mean 
field as a continuous line. For g > g c the value is 1/2, 
that is, the undecided agents have positive and negative 
a's indistinctly, but for sufficiently small g-values most 
undecided agents have a < 0. Such phase change be- 
havior is also structural as evidenced in Fig. [2^c) where 
the average cluster size ((s)) remains constant around the 
critical value g c . Here we have defined a cluster as a set of 
connected agents having opinions of the same sign, inde- 
pendent whether they are decided or not. The maximum 
observed in Fig. [2^c) and likely the minimum in Fig. [2^a) 
are due to relaxation problems. This can be understood 
due to the dynamics being stopped at a fixed time for all 
g, which is not enough to reach the asymptotic state. 
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FIG. 3: (Color online) Left panel: Time evolution of the abso- 
lute value of opinion averaged over 100 realizations of the ring 
of N = 5000 agents for g — > oo. The corresponding mean field 
prediction is shown as a continuous line. Right panels: Time 
history for a sample of 40 agents in a single realization for 
three different values of g — 7 x 10 2 , 4 x 10 3 , 4 x 10 5 (ordered 
from top to bottom), chosen to correspond to g < g c , g ~ g c , 
and g > g c , around g c obtained from Fig. [2] The vertical 
dotted lines indicate the moments when there is a change of 
time regime in the dynamics, as explained in the text. 

On the left hand side of Fig. |3]we present the results for 
the relaxation dynamics, by plotting the average absolute 
value of the state variable function of time, when 

the exchange process is off (g — > oo). The results are 
averages over 100 realizations. It is clearly noticeable 
that there are three different time regimes, seen as an 
s-shape curve and predicted by a mean field treatment 
(see Section IV). Up to around t ps g c (in units of dt) 



the evolution of (\x\) is very slow; between this value and 
t ps 2.1 x 10 5 the curve is concave upwards; finally for 
long times the variable approaches the asymptotic value 
xum — 1 very slowly and the curve is concave downwards. 
Only the evolution up to 10 7 time steps is shown, where 
fa 0.9 is reached. 
On the right hand side of Fig. [3] we show the time his- 
tory of a sample of 40 agents for a given realization and 
different values of g. If g < g c (top plot) the exchanges 



happen in the initial time regime, up to t ps g c . There- 
fore the evolution of the system cannot reach the relaxed 
state, and frustration appears in the form of indefinitely 
undecided agents. On the other hand, if g > g c (mid- 
dle and bottom plots) the relaxation is already advanced 
when the exchanges are carried out, which contribute to 
further relaxation. The fact that there is a minimum in 
the number of undecided nodes is due to the slow relax- 
ation for larger g values. 

The separation of the three time regimes mentioned 
above is even clearer in these plots: 1) in the first time 
regime practically all the agents change their opinions 
slowly regardless of the value of g, so t ~ g c can be recog- 
nized as the characteristic time for the first agents in the 
network to get decided; 2) in the second time regime the 
dynamics speeds up exponentially and most agents get 
decided (with individual trajectories getting smoother as 
g increases due to less frequent exchanges); 3) in the 
third and final time regime only some agents remain un- 
decided, which for g < g c occupy frustrated regions in 
the network and will be undecided indefinitely, and for 
g > g c will get decided after a large but finite amount 
of time. These remarks are supported by the analytical 
treatment in Section HVI 
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FIG. 4: (Color online) In the left hand side we show a semi- 
log plot of the cluster size distribution for three values of 
g = 7 x 10 2 , 4 x 10 3 , 4 x 10 s (green squares, purple circles and 
light blue diamonds respectively). The distribution for the 
initial random ring is shown for comparison as a dotted blue 
line. The panels on the right hand side show the number dis- 
tribution of a for 100 realizations of the same three values of 
g (ordered from top to bottom) and for three different ranges 
of cluster size s: red circles for s 6 [1,5), green squares for 
s £ [5, 15), and blue diamonds for s 6 [15, 5000]. 

The cluster size distribution in the initial random ring 
goes as P(s) = 1/2 S (e.g. the probability of having s con- 
secutive agents with the same sign of their initial opin- 
ions) and it is shown as a dotted blue line in the left panel 
of Fig. HJ The cluster size distribution after 10 7 trans- 
actions (also shown) turns out to change very little with 
g, although it is quite different from the random value. 
The phase change behavior seen by using (n un d) as an 
"order parameter" is also reflected in the preferred value 
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of a for clusters of different size. In the right panels of 
Fig. [4] we show this effect by plotting the number distri- 
bution of a (N((a))) for the three values of g, i.e. g < g c , 
g ~ g Cl and g > g c , and for three different ranges of clus- 
ter size. Observe that for g < g c small clusters (of size 
1 to 4) are composed mainly of agents with negative a, 
clusters of medium size (5 to 14) present a bimodal distri- 
bution of negative and positive values, and large clusters 
(15 to 5000) have agents with positive a. For g around 
the critical value the picture changes dramatically and 
large clusters start having agents with a < 0. For g > g c 
the number distribution of a approaches a Gaussian form 
independently of the cluster size. 

So far we have considered the case where the short 
range interaction deals with nearest neighbors only (m = 
1). We have also studied the situation in which the short 
range interaction includes the second neighbors, i.e. m = 
2. In this case a phase change behavior is also clearly 
visible in the number of undecided agents as a function 
of g, though the position where it appears has moved 
slightly. 



IV. MEAN FIELD CALCULATIONS 

In this section we shall investigate the peculiar features 
of the phase change exhibited by our model, namely, the 
reasons why some agents are undecided, the factors that 
determine the average cluster size, and the peculiar dis- 
tribution of agents with negative attitude parameter a 
in the final network configuration. For this we shall per- 
form a linear analysis of the dynamics and introduce some 
mean field ideas that may help understanding the role of 
the different time scales and their effect on the structure 
of the network, in particular the role of parameter g. 



A. Linear analysis 

The quantity that here plays the role of "order param- 
eter" is the number of undecided agents. However, in 
any long but finite numerical calculation, out of the to- 
tal number of agents that appear as undecided (see the 
purple squares in Fig. [2ja)) only a fraction will remain 
undecided forever. We shall investigate first the circum- 
stances that prevent agents to reach a limit opinion. 

There arc different scenarios depending on the values 
of f s and fi in Eq. [TJ When the long range term ,/) = 0, 
Eq.[T]has a simple exponential solution and the only situ- 
ation that prevents the limit value sgn(a;i) to be reached 
is when the short range term / s < 0, so the agent remains 
undecided forever. If f s = and /; 7^ the solution is 
linear in time and the agent will eventually reach a limit 
value. Notice that for m — 1 (i.e. short range interaction 
with nearest neighbors only) this situation corresponds 
to an agent at the border between two groups of opposite 
opinion, and once that agent becomes decided the border 
is displaced by one site. Since the ring is symmetric, the 



net displacement of the border will be zero, and this will 
give a characteristic cluster size. 

Eq. [T] exhibits various fixed points, on top of the limit 
values Xi — ±1. For each agent i there is a fixed point at 
\xo,i\ < 1, where 



xo,., 



flOH 
fs 



(6) 



If both f s and fi ^ 0, one can perform a linear stability 
analysis around the fixed point of Eq. [6] Then, agent i 
is considered indefinitely undecided if this fixed point is 
stable, that is, when the real part of the eigenvalue 



Xi = 



d(d t Xi) 



dxi 



= {x)f ) [29{x Q , i )-l], (7) 



is negative. In this equation 9{xq^) is the Heaviside 
step function. It should be noted that the occurrence 
of 3?[Aj] < is extremely rare without an exchange pro- 
cess (e.g. for g — > 00). The reason is that sgn(x ,i) must 
be opposite to the sign of {x)^\ which has to be different 
from zero and eventually ±1. This means that the agent 
is embedded in a very adverse environment of immediate 
neighbors, a situation not favored by the dynamics that 
tends to diminish disagreement between the agent and 
its first neighbors. The only possibility for an agent to 
remain undecided forever is when the magnitude of <Xi is 
large enough to hamper the dynamics. However, the a 
distribution is flat and the probability for this to happen 
is of the order of 0{1/N). 

Summarizing, an agent can only be undecided in the 
limit of t — > 00 if: 

(a) fi = and f s < 0, or if 

(b) /^/^Oand^Ai] <0. 

From the value of (n un d) after 10 7 transactions (shown 
in Fig.[2ja) as purple squares), we have tested all agents 
that fullfil any of these two conditions to remain unde- 
cided forever, and plotted their numbers in the figure as 
red circles. Indeed, the asymptotic number of undecided 
agents is nonzero for g < g c and zero for g > g c . The 
latter is in agreement with our previous mean field pre- 
diction, drawn as a continuous red line in Fig. [2]ja). 

We now investigate the form of the curve for the num- 
ber of undecided agents for g < g c , which can be esti- 
mated from the initial Gaussian distribution of x. First, 
the symmetry of sign in the distribution of a implies that 
only half of the agents are likely to have a < and thus 
be undecided. Then, since Eq. [T]has an approximate so- 
lution x — a;(0)e*/ Sc before the exchange process takes 
place, at least those agents with initial |ar(0) | < x g re- 
main undecided at t = g, where x g = e~ 9 ^ 5c . Therefore, 
the number of undecided agents as a function of g can be 
calculated from the initial distribution of x as 



n U nd{g) 



N erf(e-f/(ffcv / 2)/ % /2) _ erf ( e -i/V2/^) 
2~~ erf (1/^/2) 



(8) 



G 



where erf(x) = (2/y/n) JJ 5 e~ u du is the error function, 
and the factor of 1/2 is due to the sign symmetry. The 
result of Eq. [8] is plotted in Fig. [2^a) as a purple line, 
where the value g c « 1.7 x 10 3 has been fitted with least- 
squares technique. Notice that the agreement with the 
calculation (purple squares) is considerably good. The 
theoretical estimation of the truly undecided agents (red 
circles in the figure) is more involved, since the actions 
of exchanging become important, and this will be the 
matter of further study. 



B. Mean field for 

The time evolution of the average absolute value of 
opinion in the network when there are no exchanges can 
be understood by an estimation of the characteristic time 
(r) for the whole system to reach the limit values of opin- 
ion. This is done assuming that g > r, where one can by 
use a mean field approach similar to the one described in 
our previous model [55] . Although the network topology 
there is different, the mechanisms that result in the mag- 
netization relaxation of all ferromagnetic-like problems 
are similar. The average number of undecided agents as 
a function of time is found to be 

OwCO) = N-(N + 1) coth A ~ r ' M ' ' ' 



coth 



2r 



(9) 

where the time scale r is related to the critical value g c as 
r = g c N/40, see [25J. We now follow a procedure similar 
to that of the previous subsection, without considering 
exchange processes. Since Eq. [T] has an approximate so- 
lution x = a;(0)e t/,r , only the agents with initial opinion 
|a;(0)| > xt = e~ l l T can get decided at time t, while the 
rest of the agents are still undecided. Then the average 
absolute opinion of the decided agents is 2 J x P(x)dx 

and that of the undecided agents is 2e t ' T f^* xP{x)dx. 
By integrating the distribution of initial opinions P(x) 
we get 



<!*(*)!> = i- 



erf{x t /y/2) | 
erf(l/\/2) erf(l/ v / 2) < 



"Ml- e-^/ 2 



(10) 

Eq. [10] has been fitted to the numerical results shown 
in the left panel of Fig. [3] with least-squares technique, 
giving a value of r « 2.1 x 10 5 , which in turn corresponds 
to g c — 40t/N w 1.7 x 10 3 . This is in good agreement 
with our estimate of last subsection and with the value in 
Fig.j^a) of gat the minimum in the number of undecided 
agents after a fi nite number of transactions. Moreover, 
the slope of Eq. [To] is « 2 x 10" 6 for < t < g c , 

then it drops fast at around t r and is asymptotically 
zero for f>T. This illustrates the three time regimes of 
opinion evolution discussed in Section |III| and indicated 
in Fig. [3] as vertical dotted lines. We detected that the 
fitting of Eq. [lO]is very good for short times but starts to 
deviate significantly for longer times, in a similar fash- 
ion as the approximation for the fraction of undecided 



agents in [55] • This is to be expected, since in this mean 
field approach we have not taken into account the effects 
produced by the random distribution of a. As a con- 
sequence Eq. [TO] relaxes faster to the asymptotic state 
{\x\) = 1 than the actual dynamics. 



C. Analysis of (f a ) and (s) 

From the subset of undecided agents after a long but fi- 
nite time we can also calculate the fraction of agents with 
negative a. For g > g c , the symmetry of sign in Eq.[T]and 
in the initial x and a distributions implies that a should 
be distributed evenly among all agents and a value of 
(/~) = 1/2 is predicted. This agrees very well with the 
numerical results, and is depicted as a continuous purple 
line in the right part of Fig. [2jb) . 

On the other hand, for g < g c all undecided agents have 
a < 0. This fact, although apparently logical, is puzzling, 
since it holds even for reasonably large values of g (up to 
10 3 ), but it can be explained as follows: After running 
the dynamics for a long time /; is a very small number, 
since the average overall opinion approaches zero, and the 
second term of Eq. [TJ is no longer important. Therefore, 
the only way that agent i avoids the exponential approach 
to a limit opinion and remains undecided is that it finds 
itself in an adverse environment, such that it is likely 
to be chosen for an exchange many times. Furthermore, 
the exchanges have to modify the tendency of the agent 
towards a given limit opinion constantly. Remember that 
the condition to be chosen for exchange is > qff, 
meaning that the opinion in the neighborhoods of agents 
i and j is more homogeneous after the exchange. 

We now show that only the agents with negative a can 
be in this situation after a large number of exchanges 
have taken place. Consider an exchange process between 
two undecided agents with the same sign in their attitude 
parameter. If a, is positive the dynamics of Eq. [T] makes 
it likely that agent i is surrounded by neighbors that 
share its own opinion, and since agent j is in the same 
conditions as agent i, one infers from Eq. 4 that m 0. 
Since there is always a possibility that The opinion in 
the neighborhoods of the two agents have opposite sign, 
from Eq. [5] one gets qjV > and thus agents i and j are 
very seldom chosen for a location exchange, eventually 
reaching limit values of opinion. On the other hand, if 
ai is negative it is likely that Xi and (a;)' 1 ' have opposite 
signs, and a similar effect in the neighborhood of agent j 
results inp-^ > 0. Since approximately half of the agents 
j share the same sign as the neighbors of i and viceversa, 
q>V rj and the exchange can be performed. It should be 
noted that a similar analysis holds when the agents have 
opposite signs in their attitude parameter, thus « 

q^ and the exchanges are not as common as when both 
ai and otj are negative. Therefore, a negative attitude 
parameter along with the existence of many exchange 
processes hampers the possibility of agents attaining a 
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definite decision. This result is plotted in the left part of 
Fig. [5Jb) as a horizontal purple line. 

As a corollary of this analysis, we can anticipate the 
structure of a typical configuration of opinions in the sys- 
tem after a large but finite number of transactions. For 
g < g c the undecided agents with a < have been ex- 
changed many times in a random fashion and thus form 
small groups between large clusters of decided agents, as 
confirmed by visual inspection of single realizations of 
the numerical calculations. For g > g c most undecided 
agents have / s = and evolve slowly and linearly towards 
a limit opinion, therefore they should be at the borders of 
clusters with different definite opinions. This is truly the 
case, as can be seen by comparing the purple squares on 
the right hand side of Fig. |2ja) with the corresponding 
purple squares of Fig. [2jc). Indeed, the average cluster 
size is (s) 10 for large g, and the number of undecided 
agents detected in the calculation is (n un d) ~ 500, which 
is approximately NJ (s) . 



V. DISCUSSION 

In this paper we studied the coevolution of opinions 
and the embedding of individuals in their environment. 
For the opinion dynamics we adopted earlier introduced 
continuous state variable equations 28J, that include 
short range ferromagnetic interactions for describing ho- 
mophily between neighboring agents, and long range in- 
teractions for describing how the overall mood of the ma- 
jority affects the agent modulated by its attitude parame- 
ter being either positive or negative. This opinion update 
gives rise to short time scale transaction dynamics. For 
the model geometry or connectivity between agents, we 
used ring topology instead of a more complex network 
topology, we studied earlier [351 US]- The long or slow 
time scale dynamics of environment changes was carried 
out by exchanging the locations of pairs of agents. These 
two time scales are then well separated and their relation 
serves as a control parameter. 

As the main result of our study we find that by varying 
the time-scale parameter there is a phase change in the 
number of undecided individuals, which turned out to be 
mainly driven by the environment exchange dynamics. 
In order to understand this effect the following should 
be noted. First there is competing interaction due to 
the negative a's. Second, due to the asymmetry between 
the long and short range interactions (since only the lat- 
ter are considered for exchanges), this competition does 
not lead to permanent frustration, provided that enough 
time is given for relaxation. However, if the relaxation is 
hampered by too frequent changes in the neighborhood 
as well as by not allowing enough exchanges to find the 
global optimum, frustration appears as a nonzero number 
of indefinitely undecided agents. Thus the phase change 
behavior is due to the separation of time scales and due 
to insufficient relaxation. The mean field analysis which 
we performed for the system supports the above picture. 



It should be noted that there are relevant similarities 
and differences between the ring model studied here and 
our previous network model [28 . First, the transaction 
dynamics defined by Eq. [T] is equivalent in both mod- 
els, therefore producing similar relaxation processes in 
the limit of large g that can be studied with the same 
mean field approach, as has been discussed in Section [TV] 
Second, the exchange process used in the ring model is 
systematically different from the rewiring framework of 
the network model, since it changes the opinion distri- 
bution in the system by keeping the topology constant, 
thus making both models fundamentally distinct. Fur- 
thermore, the exchange process used here has some soci- 
ological background with the concept of homophily, since 
it favors homogeneity of opinion in the ring, but it is sim- 
ple enough as to allow a deeper mean field treatment than 
the one performed for the network model. 

Even though the exchange process differs from the 
rewiring scheme treated in [2 9) , both models can produce 
small clusters of agents with a < and large clusters of 
agents with a > for appropriate values of <?, as shown 
m Fig. [4] The mean field performed in the simpler ring 
topology suggests that a low time-scale ratio is responsi- 
ble for frustration that creates small groups of undecided 
agents with negative a, while a full relaxation tends to 
destroy this structure. In the case of the network model 
the rewiring rules would then be the sources of frustra- 
tion, and an appropriate value of g could freeze agents 
with a < in positions that create an heterogeneous 
community structure in the network. It is in this sense 
that a simpler though different ring model can give in- 
sight into the complex processes ocuring in our previous 
network model. 

In this paper we have analyzed different aspects of 
opinion formation and have accordingly reached the view 
that the distinction between close relations and global 
influence is an important one. While for close friends 
homophily dominates, the attitude to the overall mood 
of the society can be both positive and negative. Due to 
the latter component competing interactions occur in the 
system. Another important aspect is that - similarly to 
physical phenomena - social interactions may take place 
at very different time scales. If the separation of time 
scales is strong, the system can relax and frustration be- 
comes irrelevant, leading to a diminishing amount of un- 
decided agents. When the two time scales approach each 
other the possibility of full relaxation vanishes, then frus- 
tration appears and a finite amount of undecided agents 
remain. We think that the phase change due to a vari- 
ation in the separation of time scales is an interesting 
effect even beyond the present context. 
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